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INTRODUCTION 


The Apollo Telescope Mount Cluster (ATM, fig* l) represents the first pro- 
posed implementation of a large space station capable of supporting man and 
experiment packages in a near-earth orhit^r extended periods of time (ref. 1). 
Primary attitude control to 0.1° and 1.0 min/sec will he provxded hy a control 
moment gyro (CMG) system with gravity gradient momentum desaturation capability. 

This paper presents a discussion of theoretical and experimental research 
on the stability of this proposed CMG attitude control system. The system 
experimental data were generated with the Langley Real-Time Digital-Hardware 
Spacecraft Simulation (ref. 2). The simulation presently \ises three prototype 
1,000 ft-lb-sec CMG’s (fig* 2) and associated control electronics, and a 
CDC 6600 digital computer programed in real time for vehicle dynamics, reared 
ATM digital control logic, vehicle sensors, filters, and vehicle control logxc. 
The vehicle modal data and design- gains used were obtained from reference 5 
and will be updated as new data become available. 

The theoretical and experimental results of this work are divided into 
three areas. First, the unit CMJ mechanical and electronic characteristics ^ 
are discussed and CMG gimbal rate or precessional output torque phase and gain 
data are presented. Second, the momentum feedback or inner loop is ^scussed 
theoretically in terms of required CMG bandwidth and experimentally in terms 
of resulting moment -to -moment output phase and gain and step-function response. 
Finally, vehicle or outer loop stability data are presented as a function of_ 
CMG gimbal orientations and inner and outer loop gains. In all cases, mo^fi- , 
cations to the flight design are employed only when required to achieve stable 
low limit cycle response. The fine-pointing mode only has been evaluated at 
this writing. The effects on experimental data of simulation computer hardware 
have been minimized but are, of course, present to a certain degree. 
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square matrix of rank three (3) relating control moment G 
to angular momentum error He 

element of matrix Qa] on ith row and ;jth column 
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absolute value of element 

CMl generated control moment 

control moment component 

(Mr gimbal angle rate transfer function 

total spin angular momentum of the three (Ml's 

constant spin angular momentTun of single (Ml 

spacecraft moment of inertia 

imaginary part of s 

inner loop compensation filter gain 

ratio of simulated spacecraft attitude error gains to design 
attitude error gains 

CMG gear train torsional stiffness 

steering law forward loop gain 

maximum value of Kgj^ • H at which inner loop is stable 

spacecraft attitude error gain for ith axis 

spacecraft attitude rate error gain for ith axis 

disturbance moment components 

commanded control moment 

filtered commanded control moment 

any integer 

real part of s 

Laplace transform variable 

time 

angular frequency of bending filter 
spacecraft reference axes 

operator, operates on filtered. commanded control moment to 
yield commanded CMG angular, momentum 
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a outer glmbal coigle 

p inner gimbal angle 

A CMG gear train deadband 

0 ^ spacecraft attitude error angle of rotation about ith 

spacecraft axis 

p bending filter parameter 

0 phase angle of CMl gimbal rate relative to commanded gimbal 

rate 

ft spacecraft angular velocity con 5 >onent 

Subscripts 

c commanded value 

e error value 

filtered value 

cyclic permutation of l,2,3i x,y,z 
limited value 

component for the x-^ y-^ or z-axis 

CMGr description number; for all gimbal angles zero, 

(Mi No. 1 H-vector points along x-axis, (Mi No. 2 H-vector 
points along y-axis, CM3 No. 3 H-vector points along 
z-axis 


F 

L 

x,y,z 

1,2,3 


Brackets Cl enclosing a square array of quantities denote a square 
matrix. Braces { \ enclosing a horizontal array of quantities (separated 
by commas) denote a row vector, md braces enclosing a vertical > array denote 
a column vector. The juxtaposition of a row vector and a column vector, in 
that order, indicates that the inner product of the two vectors is to be 
taken. A bar (") over u symbol denotes a vector. A dot (*) oyer a symbol 
designates the derivative with respect to time. 


* 
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MSFC AMD LRC COORDIMATE SYSTEMS 


Figure 5 illustrates the orientation of the CMG*s and the definitions of 
the CMj gimhal angles for the MSFC and the LRC coordinates relative to the 
vehicle axes} the axis of minimum moment of inertia (workshop-CSM axis^ 
fig. l), Zj the telescope viewing axis (LM/ATM axis, fig. 1)> and y, the axis 
normal to x and z forming a right-handed system. 

The outer gimhal angles and the inner gimhal angles of the CMG’s are 
labeled V and "p," respectively, with a numerical subscript denoting the 
relevant CMS. The 5-CMS system generates control moments by means of exchange 
of angular momentum, with the remainder of the spacecraft. 

The proposed flight system conforms to the MSFC orientation} however, much 
of the data presented in this paper pertain to the LRC-CMS orientation because 
the digital simulation program in LRC coordinates had been thoroughly checked 
out, while the program in MSFC coordinates had not. The data pertaining to^ 

MSFC coordinates presented herein were acquired after the MSFC coordinate digital 
simulation program was checked out. The control systems employing the MSFC and 
LRC coordinates are sufficiently similar that most conclusions for one system 
are valid for the other. . 


CMG UNIT CHARACTERISTICS 


» The control moment gyro (CMG) (fig. 2) is a double-gimbaled momentum 
siferage device. A 24-inch-diameter rotor spinning at 11,450 rpm provides 
1,000 ft-lb-sec of momentum. The rotor is located within the inner gimhal and 
is powered by a three-phase induction-type synchronous motor. Excitation for 
the motor is supplied by either a 400 Hz supply or a variable frequency inverter 
using 28 V dc . 

The CMS unit weighs 209 pounds without the outer frame which serves as a 
movinting device. The CMS's maximum dimensions are 39*2 inches across the outer 
frame n-nd 39 inches across the outer gimbal assembly, including its acttiator and 
sensor packages. 

Located on each gimbal are an act\iator and sensor package. The actuator 
package contains a 7 ft-lb dc brush torque motor, a dc brush tachometer, and a 
two-stage planetary gear train with a ratio of 67*6 to 1. The sensor package 
consists of a second dc brush tachometer, a linear potentiometer, and sine- 
cosine resolvers. Both gimbals are fully slipringed and have complete rotation 
capability. The associated electronics provides the CMtJ with gimbal servo loops, 
rotor drive power, and system malfunction indicators such as excessive rotor 
bearing temperatures, rotor speed variations, and electronics failures. 

Figure 4 illustrates an operational block diagram of the CM3 gimbal rate 
servo loops. Command gimbal rate signals are sensed at points Aq, and Ap, 
compared with torque motor tachometer signals and the difference or error signal 
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nulled "by the torque motor transmission gimhal control loop. This results in. ia 
pseudo-control of the actual gimibad rates at points ^ and Bg» The most 
obvious feature of the servo loop is that the rate of the gimbal is not the 
controlled element but rather the rate of the torque motor. The contractor's 
design philosophy consisted of controlling the motor's rate and assu m i n g that 
the gear train (Qi^Co) between the motor and gimbal would serve as a rigid link 
with, of course, the linear reduction in rate due to the gear ratio. Gear train 
con 5 >liance, 1/%, and deadband, A, do exist, however, and have u significant 
effect on performance, as will be illustrated later. Another (MJ feature shown 
in figure 2 illustrates the momentum coupling (Do^,Dp) between the gimbals which 
causes the gimbal rate loops to interact and is, in fact, the main gyroscopic 
effect causing the predominant control torques felt by the spacecraft. 

Performance of the servo loop in CMG hardware indicated a major difference 
in frequency response between the torque motor and the gimbal. Althou^ a 
bandpass of 7 Hz was available at the motor with the contractor's servo loop 
and later updated at LRC to a bandpass of I 5 Hz, figures 5 through I 6 indicate 
gimbal bandwidths of only 3 to 4 Hz. This degradation in performance was found 
experimentally to result from the large deadband and a low stiffness of the 
gear train. Output deadbands as large as 0.1° and output torsional compliance 
values on the order of 1/47,000 rod/ft-lb were realized. The gear train dead- 
band affects the CJC gimbal response hy determining the actual magnitude of 
minimum sinusoidal displacement capable of being passed. When the tor^e motor 
is commanded to produce a constant msignitude sinusoidal rate with varying fre- 
quency, the corresponding magnitude of displacement also varies, proportional 
to l/o). Consequently, when the frequency of the commanded rate approaches 
5.5 to 4.5 Hz, the output displacement of the torque motor becomes equal to 
the input deadband of the gear train and no movement is transmitted to the 
gimbal. The compliance of the gear train also affects the gimbal response by 
providing a lag between input and output due to a spring windup effect of the ^ 
transmission. Due to the inability of the gear train to pass higher frequencies, 
a program is now in progress to improve or replace,with an optimal unit, the 
gear train and to attempt to compensate for the gear train electronic a l ly . 

Figures I6 through 19 indicate the gimbal rate responses of the No. 3 
CMG after the teeth of the gear trains had been plated, the output bearings 
preloaded, and with approximately 200 hours of operational time. A comparison 
of figures I7 through 20 with figures 13 through I6 illustrates the lack of 
improvement of the gimbal frequency responses. Although an initial improvement 
was realized with the reworked gear trains, apparently operational time has 
caused deterioration of the gear train's capability. 


MOMENTUM FEEDBACK LOOP 


Three 2,000 ft-lb-sec CMG's similar to the Langley prototypes axe employed 
in the ATM design system. Their precessional torque output is controlled by an 
inner momentum feedback loop and combined vehicle loop integrators. The design 
goal of this inner loop is to provide a high capability for accurately posi- 
tioning the three-CMG resultant momentum vector continuously as a function of 
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time. This loop is shown in_ figure 21 in block diagram form (ref. 4 ), including 
the vehicle loop integrators but excluding CMl characteristics. A complete 
dynamic analysis of this control scheme is virtually impossible because of its 
nonlinearity and complexity. However, quasi-lineeir analysis is strai^tforward 
where gimbal angle variations are ansumed small. A sin^lified case is presented 
where all gimbal angles are assumed near zero to allow the generation of 
required CMl phase and gain for inner loop stability at this particular orien— 
tation. The analysis proceeds as follows: 

Neglecting higher order small terms, the six 

Pic = KsL^Sy 
P2c = %L^®z 
Pjc 

<^lc = 

o^c = -KsLHex 

%c ' “^SL^ey 
where Kgj^ is the steering law forward loop gain. 

The three momentum error signals (where command momentum is equal to 
initial momentum) are 

Hex = H - (H - Hx2 + ®P3^ 

H^ = H - (H - Ho^ + HPi) (2) 

He^ = H - (H - Hai + H3a) 

Substitution of equations (2) into equations (l) yields the following, 

six relations. 

Pic = " ^1^ 

®5c “ %I^(Pl ” 

P2c = - P 2 ) , 

®lc “ %L^(P2 - °-l) 

Pjc = - P5) 

02c = KsiH( 35 - o^) 

Qjhese six equations are three sets of tvo equations with two variables. 

The first set is shown in block diagram form, in figure 22(a), where the Cm 
is included as an output rate to input command rate transfer function, (G(s)). 
This diagram can be simplified to that of figure 22(b) . The stability of the 
simplified loop is expressed by the inequality 
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command gimbal rates are 


(1) 


t 




-2Ksi;i 


G(s) 

s 


^CMG — ( 90°+n560° ) 


< 1 




For the first crossover (n = O) this inequality means that 




< 


^- 90 ° • 


>1 


( 5 ) 


vhere co is the CMG 90° phase lag frequency (rad/sec) and G(o)_^o) 

the CMG gain at this frequency. An average value of (Ks]^)niax using 
equation (5) for the six initial CMG servo loops was, K3LH = 2.5. 


is 


Initial simulation runs of the ATM attitude control syst«n with three LRC 
CMG’s have been completed. Typical initial stability data for the proposed 
flight system are plotted in figure 23 as stable (*) or. unstable (□ ) points 
for a set of inner loop (momentum feedback loop) gains, Kg^H, and outer loop 
(vehicle loop) attitude gain, %, expressed as a fraction of the AOM design 
gain. Constant v^icle damping was the criteria for selecting v^icle rate 
gain. Stable inner loop gain, %i,H, was consistently an order of magnitude 
lower than the design value, KjgL® ~ 10 .5. The effect on stability of 
increased CMG servo loop gain and increased noise filter bandwidth is also 
shown in figure 23* 

The primary cause of the low inner loop gains is gear train deadband and 
complianqe. Gear plating and gear bearing preloading were accomplished and 
figure 23, a plot of average CMG gimbal rate limit cycling (fine-pointing 
mode - no disturbances acting) for increasing inner loop gain, illustrates the 
Improvement in performance. 

The resulting simulation was not considered acceptable for further 
planned maneuver, acquisition, and crew motion investigations, however, since 
limit cycle magnitudes were still too large for acceptable system life and 
acceptable vehicle pointing accuracy. A compensation-filtering scheme was 
incorporated into the simulation where the actual CMG gimbal positions used in 
calculation of system momentum sure altered on a high frequency basis, but where 
the dc signal of the compensation filter remains the actual gimbal position. 
This is accomplished by adding the CMG command rate to the gimbal position and 
filtering the sum. The signal 

^'gimbal ^command /g\ 

1 + KB 


is substituted for ®gimb8LL This eliminates high frequency limit 

cycling caused by phase lag associated with actual gimbal position, ^gimbal# 
since for high frequency 
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°’glmbal ^^oamand ^ 


1 + KS 


;°^onunGmd 


and also alloys for continuous updecte of gimbal j)oaltlon. Since for low 
frequency • • ; 


^^gimbed '** ^command 
1 + KS 


“gimbal 


(7) 


( 8 ) 


Figures 25 and 26 show CMG operation with and without the above compensa- 
tion filter as the CMG's maneuver the ATM 0.05° on all three axes. No 
noticeable inner loop limit cycles are present after the maneuver is complete. 
Further data (all analog simulated inner loop) have shown that a compensation 
gain""as low as K = 2.0 with I%lH = 20.0 is possible with adequate gain 
margin. 

Figures 27(a) through (m) illustrate CMG inner loop response to a step 
command moment of 100 ft -lb. The figures show the improved inner loop response 
as the compensation-filter gain, K, is increased. Figures 28(a) and (b) 
present CMG inner loop response to a sinusoidal command moment of 100 ft-lb 
at 1.0 cps. The figures illustrate the decreased resonance and Increased 
phase associated with the compensation filter. All data are for the Integrated 
three CMG system hardware with the initial digital inner loop. 

The simulation uses CMG gimbal potentiometers rather than gimbal resolvers 
as in the ATM flight system, however, and is not in itself directly applicable 
to actual ATM CMG system hardware. The following describes a proposed method 
for using the same logic as that Txsed in the LRC simulation but specifically 
adaptable to the ATM hardware. 


Proposed Compensation Filter 

The above logic for the H-vector system change is based on substitution 
of the function 

ii^ti _ Qgimbal ^command 
1 + KS 

into the system momentum calculations rather than the gimbal potentiometer 
reading, Op, directly. Sin "a” and cos "a" are thus smooth stable variables 
and inner loop limit cycling is eliminated. 

For the ATM hardware, gimbal sine and cosine resolvers are used and thus 
Ogimbal is not available. However, the substitution of the functions 

"sin a" = “^r “)r^command 


1 + KS 



and 


tt 


cos a*' = 


(cos g)^ - (sin 
1 + KS 


vhere (sin a)^ = gimibal sine resolver signal and (cos a)^ = gimbal cosine 
resolver signal^ is equivalent to the correction used in the LRC similation; 
using only signals available in the ATM flight system. The LRC simulation is 
presently being modified to include glmbal resolver signals rather than glmbal 
potentiometer signals to validate the above corrections experimentally. 

The equivalence of both corrections can be shown by e3q>anding the trigono- 
metric glmbal angle relations as follows; 


Assuming that 



Kac 
1 + KS 


small 


then 


fl 


Sin a 


tt 


Also assuming th^t 



slnfi5ia5ai\ ~ 'Vlj.-bal 
\1 + KS/ 1 + KS 

cos/jSim^^ ~ 

\1 + KS/ 1 + KS 
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then 


sin a” = ^‘’glm'bal ^ cos ttprjjnhal 

1 + KS 1 + KS 1 + KS 


Finally, the variation between 


°^glmhal 
1 + KS 

and cos is asstuned not to he significant in the second or correction term 
(a fact to he experimentally verified). Therefore, 


g-gimhal + °^gimhal)(^c) 
1 + KS 

(sin o.)]» + (cos c,)pKaQQ]jj]iQQ 2 id 
1 + KS 


The analysis for "cos a” is similar and will not he included. 


CMS SYSTEM OR OUTER LOOP STABTLITY 


Digital simulations of the ATM-CMO control system have demonstrated space- 
craft attitude (outer loop) instahilities at certain (MS orientations. A 
summary of digital simulations of acquisition maneuvers at unfavorable initial 
CMS (noncoUinear ) orientations and some typical data traces are presented. 

Ry means of a linearized analytical treatment of the control system, minor ^ 
proposed modifications of the system 8o:e developed to stabilize the vehicle. 

A sTinmary of digital simulations of acquisition maneuvers at unfavorable 
initial CMS orientations incorporating one proposed modification of the control 
system are included herein. 

The linear analysis of the CMS control system is developed below, incoipo- 
rating both the MSFC and the LRC CMS coordinate systems (fig. 3)* 


Analysis of the H-Vector Feedback Control System 

Figure 29 schematizes the (MB control system in LRC coordinates, with 
simplifying approximations made for small vdiicle attitude and attitude rate 
errors . The simplifying approximations and the differences between the MSFC 
orientation and the LRC orientation control systems are noted in the discussion 
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Immediately following. V^icle flexibility is neglected for the linear 
analysis but included in all simulation data. 

The Euler cuigles of spacecraft attitude error, 9x, 6y# 6z trans- 
formed into commanded moment components, Mcx# Mcy# ^y means of attitude 

gains, Kx, K^, Kg, and attitude rate gains, Kpx^ K^y^ Kpg. The commanded 
moment is filtered by a fourth-order filter designed to isolate the CMS’s from 
frequencies associated with v^lcle flexibility. Each component of the 
filtered command moment is limited to ±200 ft-lb and the negative time Integral 
of the resulting moment vector is taken, with the constant of negative integra- 
tion equal to the initial value of the total spin angular momentum vector of 
the CMG’s, yielding the commanded momentum vector, Hgy, Hgg. The 

momentum error vector, formed by subtracting the total spin momentum vector of 
the CMG's from the commanded vector, generates, via the steering law, the 
commanded inner and outer glmbal precession rates . The steering law in MSEC 
coordinates is ; 


oic = KsL(-Hei sin Oi - Hej cos Oi) 


(9) 

Pic = K 3 L(-Hei cos sin P^ + Hgj sin cti sin p^ - cos P^) 


Before the commanded gimbal rates are fed to the CMG drive system, the gimbal 
rates are limited to the range ±5 .5® per sec . The outer and inner gimbsil 
angles are limited to the ranges ±175P and ^0° or ±80O, respectively. In the 
MSEC coordinates, the limits on the outer and inner gimbal angles are il?5°, 
and ±80°, respectively. The computer simulations in MSEC or LRC coordinates 
have the appropriate limits on the gimbal angles and rates and have the gyro 
transfer functions, Gj_a,(s) and Gip(s) equal to unity. The control moments 
resulting from CMGr precession are given by: 


Gi = + dj )cos pj cos ttj + (flj + djL)sin cos Pj^ 


(flj + Pj sin aj)sin Pj + (flj^ + P^ cos a^)sin p^^ 


- (Gjj sin + ftj cos Ojj + Pk)cos Pj^j] 


( 10 ) 
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In LRC coordinates ; and 


Gi = hQ- (Gk - ^)cos Pi sin Oi - (ftj - dk)cos ajj cos p^ 

+ (ftj + Pj^ cos aj^)sin Pi - + P^ sin o^)sin pj^ 

+ (flk cos aj + flj sin aj + Pj)cos pj^ (ll) 

in MSFC coordinates . These result from 

g = -d) X H 

where g is the moment produced by one CMG,_ 15 is the precessional angular 
velocity of the CMG in inertial space, and H is the angular momentum 
associated with CM} spin. 

To simplify the analysis, the vehicle attitude rates, 
ignored in equations \l0) and (ll) with little error, as is shown by the 
following order ” 0 f -magnitude calculation. It is assumed that all vehicle 
attitude rates are initially zero. 


G ~ H(ft + d) 


where G is the order of magnitude of the control moment, H = scalar angular 
momentum of a C!MG, G ~ vehicle attitude rate, d ~ glmbal precession rate. 
Note that 



where I 


vehicle moment of 


Inertia 



slug-ft^. 


Hence 


Q, 


JS_ 

IS 


where s ~ characteristic frequency of control system ^0.5 sec"l. It follows 
that 



or 
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& 


G i±^ + Hi = ^ - G + 2000a 

10°s ® 


HG 


where 2000 ft-lb-sec Is substituted for H. 

The contribution to the control moment, G, from the vdiicle rate, fl, is 

2 X 10"5 


of the order of less than 


0.3 


G, which is much less than the total 


control moment, and thus much less than the contribution from the gimbal rates. 

In the proposed flight system, the values of the sine and cosine of the 
gimbal angles, as measured by resolvers in the CMG gimbal assemblies, are used 
in the steering law and in determining the total CMG spin angular moment^ 
components, Hx, Hy, In the idealized coii5)uter s^atioM for this 

work, the gimbtl angles are obtained by integration of the simulated gimb^ 
processional rates (relative to v^icle coordinate axes) and the sines Md 
cosines are then calculated by the computer. Neglecting v^icle attitude 
rates in equations (lO) or (ll) is equivalent to consider!^ the vehicle 
coordinate system to be an inertial frame, in which case the total spin ang^ar 
momentum of the CMG's is the negative integral of the control moment, a fact 
which is useful in analysis . 

The computer simulations utilize the Euler equation. 


Gi + Ml = IiGi + (Ik - Ij)nknj 


( 12 ) 


(where M. is the ith component of the disturbance torque) to obtain the 
vehicle aigular acceleration. For analysis of acquisition maneuvers from 
initial attitude errors much less than 1° and with no disturbance torques and 
no initial vehicle attitude rates, equation (12) can be simplified to: 


G4 = l40 


i«l 


(15) 


The fractional error Introduced by simplifying equation (12) is small and 
^ = 0 , where coupling of vehicle attitude rates into off-axis Euler 

n ^ 

ancle rates is neglected (a first order in 0 effect), n is the order of 
magnitude (~) of the vdiicle attitude rates, s ~ characteristic frequency of 
system, and 0 ~ vdaicle attitude errors . Then, 
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V 


Si - ei(o) Hi at 

= 6 i( 0 ) + -^ J J Gi at dt 

■*■1 o o 


To linearize the control system, the vdiicle attitude errors are 
restricted to a small fraction of a degree, initial vdiicle angular rates are 
assumed to be zero, and the gimbal angles, gimbal rates, and commanded moments 
ere assumed to be less them their respective limits . 

The steering law in IRC coordinates is; 

®ic = %L("^ei ““i " ®ek 

Pic = KsL(-Hei cos«Oi sin Pi + Hej cos Pi + Hek sin oi sin Pi) 

The CMG control moment components are approximately: 

Gi = H(a^ cos Pj cos aj + <xi sin ai cos Pi - Pj sin aj sin Pj 
+ Pi cos tti sin Pi - Pk cos P]j) 

in IRC coordinates. The vehicle b^avior as simulated by con^ruter has been 
observed to be insensitive to substitution of ideal CMG*s (for which 
G (s) = Gq(s) = l) for actual CMS hardware when the inner loop (with hardware) 
is stable^ (see earlier section). The assumption follows that 

Oi = die 

Pi = Pic 

The CMG control moment components become, in IRC coordinates: 

I cos^j cos Pj + sin^ai cos Pi + sin^Pj sin^j 
sin aj cos aj cos Pj 

sin Pj cos Pj sin aj - sin^Pl sin Oi cos oti 


+ cos^l sin^Pl + cos^Pjj 

- sin aj cos aj sin^Pj - sin Pi cos Pi cos <xi + sin Pk cos Pk sin aj^ 

- sin Pit cos Pit cos ttk + sin oi cos ai cos Pi 
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s imil arly j in MSPC coordinates^ 


r 


Gi = -KsIjH Hgj, 


sin^ cos + cos^ajf cos P]c 

p 

sin cos Oi cos - sin cos sin Pi 

sin ojj cos cos Pj^ + sin Pi cos Pi cos ai 

2 . 


+ cos^tti sin^Pi + sin^ajj sin^Pjj + cos'^^Pj 

- sin Pjj cos Pjj sin Oj^ + cos aj sin Pj cos Pj 

- sin Ojj cos oj^ sin Pj^ - sin aj sin Pj cos Pj 

Expanding equations (l6) and (iT), define the matrix [V] by; 




^11 

Ai2 

^15 

= ^Hex^ 

Hey, Hez^ 

A21 

A22 

A25 



> 

A52 

^5. 

“ ^^ex> 

Hey# QQ 




(IT) 


(18) 


A 

for both MSEC and LRC coordinates . A 21st order characteristic equation of the 
control system ceui be derived with no additional simplifying assumptions • 
Neglect of the components of the matrix [Aj] off the principal diagonal results 
in the seventh order characteristic equation derived in the following, using 
figure 29 . The predictions of vehicle behavior based on this simplified theory 
are then verified by observed computer simulated vehicle behavior. Heye, 

( 19 ) 


and 

= (ao) 


by d^inltion. Note that 


1%1 
S-. 


(21)<}. 

15 


I 


and 


-Hl=% 


( 22 ) 


hence 


or 


Gi = Aii 




=FTT- 


li 


For KslH> 0 and \^^\, |pj|, |Pk|<|^ 


Aii <0 


and 


G. = 


s + 


l^iil 
l*iil 


Mpi 


The output of the bending filter Is 


(25) 


(24) 


%1 = Mci 


s**" (4p| + 2) ^ 4p_ 

—— + s 5 + + s + 1 

wj; wj wf ^1 




(25) 


•where 


^cl “ "%®1 " %!®®1 


(26) 


Therefore 


0i = 


I All 


s + 


|Alll 


{-Kl01 - %is9l}r^ 


h.M1£1s ti£is + i 
wf '^l 


( 27 ) 
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Substituting the expression 


Gj = 

Into (27) and arranging In standard form yields 


g7 + 

"|Aiil h , ‘‘PiiJ 

s6 + 

(4pf + 2)lj ^ 4pjli 






( 28 ) 


4pili (4pf + 2)li 

I ■ 




Wf 


s^ + (Ii)s5 + (|Aii|li)s^ 


+ (%llAiil)s + %|Aii| =0 (29) 

The above Is valid for either MSEC or LRC coordinates. 

Equation (29) vas solved numerically, varying the parameter. |Aii| , for 
present design values of the functions 1^, pj^, and vhlch are: 

Ix = 0.4555 X 10^ slug-ft^ 
ly = 5.475 X 10^ slug-ft^ 

I, = 5.235 X 10 ^ slug-ft 2 

= 4.26 rad/sec 
Wg = 1.5 rad/sec 
Wj = 2.6 rad/sec 

Pi = 0.5 
P2 ® 0.5 
Pj = 0.5 
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= 1.475 X 105 ft-lb 

Ky = 2.95 X, 10^ ft-lb 

Kz = 7075 X 10^ ft-lb 


Kpjj = 2.95 X 10^ ft-lb-sec 

Kt>, = 1.1 X 10^ ft-lb-sec 

By 

%z = 2.46 X 10^ ft-lb-sec 


The parameter |A^i| was varied from the maximum possible value when the 
steering law constant K^t is the design value, 0.005^ slug“^-ft“2^ to • 
values which resulted in v^icle instabilities (positive real parts of any 
solution for s). In LRC coordinates, the element is, from equation (l6), 

Aii = -Kgj^H(cos2aj cos Pj + sin2ttj^ cos + sin2pj sin^j 

+ cos2a^ sin2pj|^ + cos2pj^) (jO) 


The limits on the inner gimbal angles, P^, are plus and minus an acut^ 
angle, and so, cos Pi > 0, cos Pj > 0, and cos p^ > 0. Thus all the terns 
in parentheses in equation (50) are positive, and the minimum value of jAiij 
results when ai = 0 or it, aj = Pi = 0, Pj = 0, ±Pk = PL, yielding 

lAiil equal to KslH cos2p^ sec-1. In MSFC coordinates, the element A^t is, 
from equation (17) • 

Ail = -KsLH(sin^o^ cos Pj^ + cos^o^ cos Pj^ + cos^a^^ sin^Pt 

+ sin2ttj^ sin2pjj + cos2pj) (51) 


The TnlnlTmiTn value of [Aiil results when ^i = 0 or it. 

Pi = 0, Pk = 0, Pj = ±Pl> yielding |Aii| equal to KslH cos2pL sec"^. 

For Pl = *80 , I Aiil a minimum value of O.517 sec in either MSFC or 

LRC coordinates. For Pl - ^ 0°, I Aiil has a minimum value of 1.25 sec“l. 


Figures 30, 31, and 52 display the real parts of the solutions of 
equation (29) for s for 0x> 0y# and 0z> respectively. The least damped 

root, for all axes, rapidly approaches zero damping as |Aii| drops. from 
about 5 sec“^, indicating deteriorating vehicle control capability. Figure 35 
is a plot of the least damped root for all axes in the neighborhood of neutral 
stability (zero damping). The linear model predicts neutral stability at 
| An I = 0.963 sec“^ or | App | = 0.825 sec"^. or IA55I = 1.10 sec“^. Digital 
simulations of acquisition maneuvers from initial spacecraft attitude errors 
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of 10"^ rad (~ 0.006°) on all axes (vhich is in the linear range) with ideal 
gyros (go,(s) = SpCs) = l) demonstrated neutral stability at the following 
initial glmbal angles, in LRC coordinates: 

0x neutral stability; Og = Pj =68,8°, other gimbal angles zero 

0y neutral stability; Pi = 71*0°, other gimbal angles zero 

02 neutral stability; = 90°, P2 “ 6o.2°, other gimbal angles zero 

where Pl “ The linear model predicts neutral stability at the following 

gimbal angles: 

0x neutral stability; ~ 90°, Pj = 72.^°, other gimbal angles zero 

0y neutral stability; = 90°, Pi = 75*7°, other gimbal angles zero 

0^ neutral stability; = 90°, 02' = 71*2°, other gimbal angles zero 

by virtue of the fact that | An I = Kh]^ cos^(72.4°) = 0.963 sec"^, 

IA22I = cos 2(73.7°) = 0.825 sec-l, [Ajjj = Ksi^I cos2(7l.2°) =1.10 sec“l, 

respectively, where Kjgj^ is O.OO525 X 2000 = 10.5 sec"^, the design value, 
figures 34(a), 34(b), and 34(c) show computer simulated vehicle attitude for 
three acquisition maneuvers in the neighborhood of neutral stability for 

There are (MJ gimbal orientations in which the off-diagonal components of 
the matrix Q0 vanish. Some examples are, in both LRC and MSFC coordinates: 

(a) all gimbal angles =0 

(b) . = ^90°> other gimbal angles zero 

(c) 02 = *90° > other gimbal angles zero 

(d) 05 = ±90°, other gimbal angles zero 

as can be seen from equations (16) and (17)* • As a check of the linear model of 

the control system vinder conditions where the complicating effects of the off- 
diagonal elements of the matrix |aT are not present, computer simulations of 
acquisition maneuvers were run ar gimbal orientations that diagonalize lAj. 

The value of Kg^H was varied in order to induce neutral stability and 
instability. The following summarizes the results. 
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NEUTRAL STABILITY, (LRC COORDINATES) 


Initial 

glmbal 

angles 

Initial nonzero 
attitude angles 

KslH 

simulated 

KslH 

predicted by 
lineeir model 

02 = 90°, 

others zero 

0 jj = 10 ”^ rad 

■ 

0.9625 

02 = 90°, 

others zero 

©X = 0 y = 0 z “ 

■ 

0.825 

ai = 90°, 

others zero 

— k 

0 g = 10 rad 

■ 

1.1 


Figures 55(a), 35(b), and 35(c) are recordings of the simulated vehicle attitude 
for these acquisition maneuvers. The predicted value of for neutral 

stability is fo\ind by means of equation (30) and figure 35 as follows: 


For . 

tti = 90° 

and others 

zero. 

Aii = -3KslH 

For 

aj = 90° 

and others 

zero. 

Aii = -KslH 

For 

ttk = 90° 

and others 

zero. 

All = -2KslH 


For neutral stability in 0 x, A^ = -O.9625. 

For neutral stability in 6y, A22 = - 0.825 

For neutral stability in Sg, A^j = - 1.1 


Therefore, predicted neutreO. 

0x °2 " 90° 

0 y at aj = 90^ 
0 g at = 90° 


stability occvirs for 


and others zero when Kgj^ = O.9625 
and others zero when = 0.825- 
and others zero when KgLH = 1..1 
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Figures 56, 57, and 58 display the imaginary parts (frequencies) of the 
complex solutions of equation (29) . The negative imaglnarles of the complex 
conjugates to the roots with positive imaglnarles are not plotted. The frequencies 
of the least damped solutions are observable in computer simulations of the CMG 
control system near the condition of neutral stability with respect to the axis 
to which the solution pertains; the wave forms of the spacecraft attitude then 
persist long enough for the wavelengths on the recording chart to be measured. 
Figure 39 is a plot of the solution frequencies in the neighborhood of neutral 
stability and of frequencies observed in comx>uter simulations of the control 
system in LRC coordinates with ideal CMG's. The "coupled x" data points 
correspond to the condition 02 = P5 “ other gimbal angles zero, and 

KotH has the value required to make I A n \ equal to the value on the abscissa. 
Similarly, "coupled y" means aj = 90°, Pi = 70°, and "coupled z" means 
Op = 90 °, Pt = 70°. "Uncoupled x, y, z^ correspond to the condition 
Op = 90°, ttT = 90°, and = 90°, respectively, withj^KsjjH variable. The 
initial vehicle attitude errors are 0x = 9y = 0z = 10 rad except for the 
"uncoupled x" and "uncoupled z" cases, where the nonzero initial attitude 
errors are 0x = 10"^ rad and ©z = 10"^ rad, respectively. The te^, 

"coupled," Indicates that the off-diagonal elements of the matrix [Aj are not 
all zero. 

Table II is a summary of simulated acquisition maneuvers in MSFC 
coordinates, incorporating the isogonal correction logic. The isogonal scheme 
tends to minimize the inner gimbal angles and to space the individual CMG spin 
angular momentum vectors symmetrically with respect to the total CMG spin 
angular momentum vector. The simulations of table II are outside of the 
linear range of the control system; the filtered commanded moment limit and 
the gimbal rate limit are both invoked. 

Equation (24) and the linear model Indicate increasing lag of the control 
moment, Gi, with respect to the filtered commanded moment, as the 

magnitude of the matrix element An decreases. One way to tiy to avoid the 
instability associated with low |Aii| is to raise the lower limit on |Aii| 
by increasing i KslH cuid/or narrowing the limits on the inner gimbal angles . 

The upper limit on Ksj^H is set by the requirement of inner loop stability. 

The inner loop is sketched in figure 29 and discussed in the previous -section. 

The penalty paid for narrowing the inner gimbal travel is that the inner 
gimbals will hit their stops more often than before, thus reducing CMG 
redundancy and torque linearity and/or more power will be consumed in moving 
the inner gimbal away from the stops. These possible remedies will be investi- 
gated in the future. Another stabilizing modification of the. control system 
is to replace the negative inte^ation of the limited commanded moment, Mj-, by 

an operator Z, such that The following is a derivation of the 

operator, Z, utilizing figure 29, for both MSFC and LRC coordinates. With 


“ All- 


(52) 
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Mpi 

and replacing the relation, by in figure 29, 

Hei = Hci - Hi = Mpi • Z - Hi 
Gi 

= %i • Z + -F . 


Hy equations (52) and (18), 

^ = %i • Z + ^(AiiHei + AijHgj + Au^He^) 

so that 

Z = + ^ %k\ 

" Ail s^ AjJ %i ^kk %i/ 


( 53 ) 


the constant of integration being the initial value of Hi* The control moment 
becomes equal to the control moment of the unmodified control system at 
infinite I All I, that is, there is no lag from ^^i to Gi. Note that 
division by the An's is allowed because they are never zero. The control 
moment Gi is, by equations (18) and (32): 

Gi = AiiHei + AijHej + Aij^Hgk 

Figure kO schematizes the assumed linear circuit equivalent to the circuit 
of the control system employing the operator Z. The equivalent circuit has no 
inner loop, however, the idealized computer simulation represented by figure kO 
and used to test the utility of the Z operator is equivalent to the idealized 
computer simulation of the GMJ control system incorporating the inner loop and 
the operator Z, when all simulated CMG’s are fully functional. 

The real parts of the least damped roots of the sixth order characteristic 
equation for the linearized control system represented by figure 40 , neglecting 
the nondiagonal elements of the matrix A, are: 
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For Q^, -0.552 sec“^, independent of 
For 0y, -0.167 sec”^, independent of 
For 02 > -0.590 sec”^, independent of 


Comparison of these root parts vlth the real parts of the roots of the 
characteristic equation for the linearized control systan of figure 29 
indicates that the Z operator increases vehicle stability at low values 
|Aii| . Idealized computer simulations demonstrate that the Z operator 
improves the CMS systan control capability overall. Table III is a summary 
and comparison of idealized digital simulations of v^icle acquisition 
maneuvers performed by the CMG control system with and without the Z. operator. 
The simulations, in LRC coordinates, are outside of the range of the 

control system because gimbal rates and filtered commanded moments are confined 
by the limiters. 
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CONCLUSIONS AND RECOMMENDATIONS 


The conclusions and recommendations for this paper are divided into the 
three ATM-CMG systoa areas; (l) the CMS servo loops, ( 2 ) the ATM H-vector 
inner loops, and (5) the controlled ATM vehicle attitude loops. 

1 . The main CMG problem areas consist of developing and inclement ing a 
responsive and accurate gimbal rate servo loop consisting of optimized 
mechanical and electronic components. Mechanically, optimization consists of 
the utilization of the most efficient and linear components possible such as 
brushless dc torquers, low- compliance — low-deadband gear trains, etc. 
Electronically, optimization consists of designing servo loop con5>ensation 
which will allow wide bandwidth control of the gimbal rate itself with the 
ability to cancel nonlinear mechanical effects introduced by elements such as 
gear trains. As demonstrated by the performance of the CMG's in spacecraft 
simulations, the ability of CMl*s to accurately control the attitude of a 
spacecraift is directly dependent on the responsiveness and accuracy of the 
gimbal rate loops. Consequently, it is highly desirable to develop servo 
loops approaching the ideal time optimal response characteristics. A CMl unit 
optimization program is underway at Langley Research Center which will provide 
time (and energy) optimal second and third generation CMG hardware for attitude 
control systems of this type. 

2 . The concept of H-vector or CMG momentum feedback has been found to 
be an extremely difficult concept to implement with state-ofrthe-eirt CMG 
hardware. The LRC digital simulation has shown that for the required ATM 
design gains the proposed system will not function satisfactorily. The inner 
loop coD5>ensat ion-filter scheme discussed in this report will eliminate the 
inner loop instabilities for the LRC simulation and should have a similar 
effect for the fll^t system, but is in no way a substitute for compatible 
control logic -CMG hardware. 
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3. Analysis of a lin^r model of the ATM CMG control system reveals some, 
but not all, of the mechanisms for vehicle (outer loop) instabilities. One of 
the causes of vehicle instabilities is excessive phase lag of the control 
moment with respect to the commanded moment. The excessive phase lag results 
from low magnitudes of the diagonal elements of the momentum-error-to-control- 
moment transfer matrix at certain CMG orientations . The incorporation of a 
filter, termed the Z operator, into the control system compensates for the 
phase lag and eliminates instabilities in both the linear and nonlinear regions . 
The modified H-vector control system (Z operator modification) shows sufficient 
promise to warrant further study although the additional computational require- 
ment coupled with the present computational requirement of this system may 
exceed a rational limit. The possibility of stabilizing the spacecraft by- 
increasing KslH and/or narrowing the limits on the inner gimbal angles, as 
an alternative to the z operator, will be investigated in the future. Also, 
the response of the system with a CMG failed has been shown to cause outer 
loop instabilities. This is being investigated and data will be published in 
the future. 
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V. 


TABLE I CMS Constants 

i ■ ■ 

i 

! 

K » 184' 

Kl = 6.8 
L =22 volts 
K]» = 0.l6 ft-lb/volt 
Ij » 0.005 SLUG ft2 
G » 6j.6 

Kq = 47,000 ft-Ib/rad 
Kh = 0.427 

= 4.9 SLUG ft2 
Ip =4.0 SLUG ft^ 

Tf = 0.15 ft-l'b 
A = 1/573 rad 
H = 1000 ft-lb/sec 
Kg = 1.53 volts/rad/sec 
=0.1 sec 
T2 = 1*0 sec 
T3 = 0. 5 sec 
t4 = 0.003 sec 
T3 =0.02 sec 

Tg = 0.0033 sec 


TABLE II 

SUMMARY OF ACQUISITION MANEUVERS, H-VECTOR WITH ISOGONAL CORRECTION, MSEC COORDINATES 


/O - +80° 

Inlxial non-zero 
Gimbal angles^ deg. 

= goO,yS 2 = 80° 

i 0(3 = 90°,/? 2 = 75° 

1 

Initial Spacecraft 
Attitude errors, rad. 
= 9.y = rz 

.8725 X 10 "^ 

.8725 X 10 “^ 

.8725 X 10"3 
.8725 X 10"^ 

Stability 

temporarily unstable x 
temporarily unstable x 

Acquisition 1 

Time, .990, 

• ■ - — -i 

90 . 1 

75 . ■! 

85 1 

’ 60 ■■ j 

1 = 9O°,>0 3 = 75° 

temporarily unstable y 
stable 

o< 2 = 90°,^1 = 75° 

.8725 X 10 "^ 

temporarily unstable z 

65 • i 

2 ~ ^3 ~ 

.8725 X 10"3 

unstable x 

I- 

1 


.8725 X 10 ’^ 

stable 

45 i 

= 70° 

.8725 X 10”9 

iinstable y 

\ 

/3 1=/32= 

.8725 X 10“3 

temporarily unstable z 

220 ■ • j 

1 o( 3 ^ = 90° 

.8725 X 10"9 

stable 

55 • j 


.8725 X 10 "^ 

temporarily unstable y 

280 i 

o( 1 = - 90 ° 

.8725 X 10"^ 

unstable y 

• ! 

•— I 

: <^ 2 =*° 

.8725 X 10-3 

stable 

' 1 

c< 2 = 90 ° 

.8725 X 10-2 

stable 

200 • . s 

I 

ol = 90° 

.8725 X 10“2 

1 stable 

160 ■ ‘ 

— - 

.3723 X 10-2 

unstable x 


. 4 . 3 ..= r9Q — 

2 = “90° 

.3725 X 10-2 

unstable y Sc z • 







f 


TABLE III 




SUMMARY OF ACQUISITION MANEUVERS, LRC 

COORDINATES 




^ L = ± 70° 

Initial Spacecraft 

Stability 

Acquisition 

Control Law 


Initial Non-zero 

Attitude Errors, rad. 

Time, sec. 



Gimbal Angles, deg. 

9x ~ ~ 



New (as per Fig. 

41) 

c\ 1 — 

.8725 X 10-3 

stable 

30 

Old (as per Fig. 

2-?} 

o/i=.x'2=c(3 = 45° 

*8725 X 10"^ 

stable 

40 

New 


ol 2 - ^'3 = 

.8725 X 10"^ 

stable 

40 

Old 


C< 2 = 90°, (3^ = 70° 

.8725 X 10 "^ 

stable 

25 

New 


r <3 = 900,/?! = 70° 

.8725 X 10-3 

stablo, initial x overshoot 

60 

Old 


Cv’ 3 = 90°, = 70° 

.8725 X 10“3 

rsiotablo x and y 

— 

, New 


c^.-^ = 90°, /^2 - 

.8725 X 10"3 

stable, initial 2 overshoot 

70 ' 

Old 


= 90 °, f^2 = 70° 

.8725 X 10-3 

temporarily unstable 2 and y 

85 , 

New 


02 = ^3 = 

.8725 X 10-3 

stable 

. 40 

Old 


^ 2 = ^3 = 70° 

.8725 X 10 “^ 

.unstable x 

— 

! New 


^2=P^3= '70° 

-.8725 X 10 “^ 

stable 

35 

Old 


/-■ 2 = (h = '70° 

-.8725 X 10-3 

unstable x 

— 

New 


Pi = ^ 3 = 70° 

.8725 X 10-3 

stable 

50 

: Old 


Pi =.^3= 70° 

.0725 X 10-3 

unstable y 


New 


P 1=^2 =70° 

.3725 X 10"3 

stable 

45 

Old 


pi =P2 = 70° 

.3725 X 10”^ 

temporarily unstable 3 

68 . 
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I 

j 

I TABLE III, -CONCLUDED 

V . ‘ ■ I, 

■ ’ •• • 

! 


i i 

1 Control Law 

i * ■ 

i 

^h = ± 70° 
Initial Non-zero 
Gimbsd Angles, deg. 

Initial Spacecraift 
Attitude Errors, rad. 

=&y = 

Stability 

i 

Acquisition j 
Time, sec. 

1 

! , New 

4l =^2 - *7°° 

-.8725 X 10“^ 

stable 

35 

Old 

^1 =^2 = 70° 

-.8725 X 10“^ 

temporarily unstable z 

70 

New 

c/ 1 = 90°, = 70° 

10“^ 

stable 

70 

1 Old 

90°, f^2 - 70° 

10“^ 

stable, initial 3 , y, z 
overshoot 

45 

1 ' Old 

90°, ^2 = 70° 

lO”^ 

stable 

110 

Old 

o(l = 90°, (^2 = 70° 

1.5 X 10"5 

temporarily unstable 

160 

New 

^1 = 90° 

.8725 X 10“^ • 

stable 

no 

. Old 

c^i = 90° 

.8725 X 10“3 

unstable y 
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TABLE IV*-“ Conversion Table 


U. S. Customary Units 

S. I. Units 

inch 

.025^ meter 

foot-poiind 

1 . 356 Neuton-meter 

pound 

Newton 

2 

slug-foot 

2 

1.356 kilogram-meter 
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Figure 5 .-Low amplitude frequency response 
- inner gimbal CI-® #1 • 
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High amplitude frequency response 
- inner gimbal CMG #2 
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High amplitude frequency response 
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Figure 13. 


Low amplitude frequency response 
- inner gimtal CMGjjr3 
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High amplitude frequency response 
- inner gimbal CKG ?r3 
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L_,,_ ,c Low amplitude frequency response 
Figure 15. _ cKG #3 
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Figure 1(. _ gimbal CMG 


1 

FREQUENCY (Hz) 


-240 

10 


X 



>i 

I 

i 151- 

( • 

Inner Glnibal - #3 

8* pp/sec cdmaand 
o” Oiinibal Angle 
Plated Gear Train 











i 


i 

10 - 


Outer Glidbal 080 #3 
2* pp/sec cdmnapd 
p* aimibal Angle 
Plated -Geax: Train 



-210 



Low amplitude frequency response 
- inner gimbal - CM3 #3 


i 

FPJi:Q,UENCY (Hz) 


10-240 


Figure 19. 


1 


15 


Outer Gimbal cJfG #3 
8® pp/sec command 
0“ Global Angle 
Plated Gear Train 



!0 

! PHASE 
i(deg) 
;-3o 


-6o 


.-90 

^-120 


-150 


.-l80 


-210 


-20 


1 


Figure 20. 


High amplitude frequency response 
- inner gimbal O'®} /■•3 
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- COMPENSATED CONTROL SYSTEM RESPONSE 
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Figure 2.S.-UNC0MPENSATEP CONTROL SYSTEM RESPONSE 
MEASURED GIMBAL RATES FOR THREE- AX IS .05° SPACECRAFT MANEUVER 
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Figure 29.- H-vector feedback control system. 
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Figure 34.- Simulated acquisition maneuvers. 
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Figure 35.- Simulated acquisition maneuvers 
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Figure 40.- Equivalent circuit for H-vector feedback system with Z operator. 
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